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FOREWORD 


This  is  Volume  II  of  three  separately  bound  volumes  in  which  are  reported  the 
research  completed  under  the  general  terms  of  the  Office  of  Civil  Defense  Subtask 
Number  4113E,  "Sensitivity  Analysis  of  Civil  Defense  Systems  and  Components." 

The  author  acknowledges  the  valuable  assistance  of  Mr.  Philip  McMullan  in 
completing  the  research  and  of  Mr.  McMullan  and  Mr.  Robert  Brooks  in  preparation 
of  the  final  report. 


ABSTRACT 


In  order  to  perform  a  sensitivity  analysis  of  selected  parameters  of  interest 
in  Civil  Defense  systems  analysis,  probable  casualties  are  estimated  for  8  SMSA's 
over  a  range  of  fallout  environments  and  shelter  utilization  patterns.  The  selected 
parameters  are:  SMSA  population,  population  density,  and  ratio  of  shelter  spaces  to 
population;  fallout  arrival  time  and  reference  intensity;  and  restrictions  on  move¬ 
ment  of  people  to  shelter,  leading  to  varying  patterns  of  shelter  utilization.  The 
SMSA's  are  selected  by  "judgment  sampling"  and  range  in  population  from  74,000  to 
408,000.  The  fallout  environments  used  range  from  a  reference  intensity  of  600  r/hr 
and  7  hours  time  of  arrival  to  a  reference  intensity  of  30,000  r/hr  and  1  hour  time 
of  arrival.  The  movement-to-shel ter  restrictions  are:  (1)  movement  restricted  to 
the  Standard  Location  (SL)  of  residence,  (2)  movement  restricted  to  within  two  miles 
of  the  SL  of  residence,  and  (3)  unrestricted  movement  to  shelter  anywhere  within  the 
SMSA.  Also,  (4)  the  transportation  algorithm  is  used  to  determine  the  optimal 
(minimum  casualty)  allocation  of  people  to  shelter  for  each  time  of  arrival  and 
reference  intensity  combination.  This  allocation  serves  as  a  benchmark  of  ideality 
against  which  to  measure  other  patterns  of  shelter  utilization.  Casualties  are 
computed  for  each  of  the  four  movement  patterns  over  the  range  of  attack  environ¬ 
ments.  It  is  concluded  that  efficient  shelter  utilization  is  very  important  when 
shelter  spaces  exceed  population;  inefficient  utilization  can  lead  to  a  large  number 
of  avoidable  casualties.  Movement  through  fallout  for  a  short  time  in  order  to 
reach  a  better  shelter  theoretically  is  warranted  in  many  circumstances. 
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A  Sensitivity  Analysis  of  Selected  Parameters 
Based  on  8  SMSA's 


I .  INTRODUCTION 

The  purpose  of  this  volume  is  to  analyze  the  sensitivity  of  estimates  of  the 
effectiveness  of  fallout  shelter  systems  to  changes  in  movement  assumption?  and 
resultant  shelter  utilization  under  various  fallout  environments,  shelter/population 
ratios,  and  population  densities. 

Eight  Standard  Metropolitan  Statistical  Areas  (SMSA's)  were  chosen  representing 
a  range  of  population  densities  and  shelter/population  ratios.  Four  modes  of  alloca¬ 
ting  population  to  shelters,  based  or,  travel  restrictions  in  time  or  distances,  and 
six  types  of  fallout  environment  were  assumed  for  each  of  the  SMSA's.  Expected 
casualties  were  computed  for  each  combination  of  allocation  mode  and  fallout  environ¬ 
ment. 

The  results  were  then  analyzed  to  answer  these  and  other  questions: 

(1)  How  effective  are  various  pat* ?rns  of  shelter  utilization  in  each  of  the 
SMSA's? 

(2)  What  trade-offs  exist  between  protection  factors  and  movement  in  fallout? 

II.  DEFINING  THE  MODEL 

The  four  types  of  movement  restrictions  considered  in  the  analysis  are  as 
follows: 

(1)  Movement  restricted  to  within  each  standard  location  (SL)  and  no  movement 
in  fallout  necessary  to  reach  shelter. 

(2)  Movement  restricted  to  within  an  SL  or  to  SL's  not  more  than  two  miles 
(distances  are  rectangular  distances--rather  than  straight  line--measured 
between  the  centers  of  SL's);  no  restriction  on  speed  of  movement. 

(3)  Movement  unrestricted  within  the  SMSA  and  no  restriction  on  speed  of  move¬ 
ment. 
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(4)  Movement  unrestricted  within  the  SMSA  and  a  constant  movement  speed  of  two 
miles  per  hour  beginning  with  weapon  detonation;  movement  in  fallout  per* 
mitted;  casualties  minimized  using  the  transportation  algorithm  (see  below). 

Assumptions  for  all  cases  were  as  follows: 

(1)  All  standard  locations  in  an  SMSA  have  the  same  fallout  environment; 

(2)  The  fallout  intensity  increases  linearly  during  the  time  fallout  particles 
are  being  deposited; 

(3)  Maximum  Equivalent  Residual  Dose  (EPD)  is  calculated  by  a  finite  differencing 
method  which  examines  dose  accumulation  and  recovery  from  time  of  arrival 
through  the  initial  shelter  period  (see  Appendix  B) ; 

(4)  Midpoint  values  of  each  Protection  Factor  (PF)  category  are  used  in  casualty 
calculations; 

(5)  A  person  experiencing  a  Maximum  ERD  of  200r  is  considered  a  casualty. 

Calculations  of  casualties  for  the  case  with  no  movement  outside  the  standard 
location  were  straightforward.  NFSS  Phase  1  data  [Reference  1]  were  used  for  each 
SL  to  determine  the  population  and  shelter  spaces  of  each  PF  category,  including 
Category  1.  The  population  was  then  allocated  manually  on  the  basis  of  utilizing 
the  highest  available  PF  shelter  space.  Any  population  in  excess  of  the  number  of 
shelter  spaces  available  was  considered  to  have  a  PF  of  2. 

In  the  third  type  of  movement  restriction  (movement  anywhere  within  the  SMSA 
and  no  speed  restriction),  population  was  allocated  to  the  highest  available  PF 
shelter  spaces,  regardless  of  their  initial  locations  within  an  SMSA. 

In  the  case  of  movement  restrictions  of  two  miles,  it  was  necessary  to  compute 
the  distances  between  centers  of  SL's.  The  rectangular  distance  was  used  since  it 
is  reasonable  that  the  route  to  shelters  in  most  cities  must  follow  some  right  angle 
street  pattern  rather  than  a  straight  line.  The  allocation  of  population  to  shelters 
was  then  done  manually,  always  using  all  the  space  with  the  highest  PF  within  the 
two-mile  movement  restriction. 

The  fourth  case  of  movement  restriction,  limiting  movement  speed  to  two  miles 
per  hour,  was  much  more  complex.  It  was  necessary  to  include  the  possibility  of 
movement  in  fallout  in  order  to  reach  a  more  effective  shelter.  In  order  to  cal¬ 
culate  the  movement  pattern  leading  to  the  minimum  number  of  casualties,  it  was 
first  necessary  to  calculate  the  ERD  accumulated  in  moving  from  one  standard  loca¬ 
tion  (S^)  to  shelter  in  another  standard  location  (S j )  of  protection  factor  (It). 


When  the  travel  tine  between  standard  locations  was  great  enough  to  require  movement 
in  fallout,  it  was  necessary  to  determine  this  exposure  time  and  to  calculate  the 
dose  sustained  outside,  as  well  as  that  sustained  while  in  shelter. 

The  process  of  determining  the  EPD  for  movement  between  standard  locations  under 
various  fallout  environments  was  programed  for  the  CDC  3600  computer.  The  generalized 
matrix  that  was  used  for  these  ERD  "costs'*  is  shown  below  in  Table  I. 


TABLE  I 

Generalized  Matrix  for  ERD  "Costs" 


S1 

S2 

Destinations 

sj 

Origins 

PF1 

pf2 

PFl 

pf2 

•  •  •  •  • »  PF^ 

•  •  •  • 

PFk 

S1 

Slll 

S1 1 2 

S121 

S1 22 

..  ..  ..  Sj 

•  •  •  • 

siJk 

S2 

•  • 

S211 

•  • 

S212 

•  • 

S221 

S222 

.  S2jl 

•  •  •  • 

S2jk 

•  • 

•  • 

•  • 

(e. g. ,  S^ 22 
destination  2 

the"ERD"  cost  to  move  from  origin  1 

and  enter  a  PF  category  2  space) 

to 

Si 

Sill 

Sil2 

Si21 

Si22 

.  Sijl 

•  ♦  •  • 

S.  .. 
ijk 

With  such  a  matrix,  and  with  the  assumption  that  casualties  occur  when 
ERD  >  200r,  the  transportation  problem  method  of  linear  programming  was  used  to  mini¬ 
mize  the  number  of  casualties  within  each  SMSA  for  each  fallout  environment.  The 
mathematical  description  of  the  transportation  model  is  given  in  Appendix  B. 


III.  DATA  FOR  EVALUATION 


A.  Selection  of  Local  Areas  (SMSA's) 


Eight  SMSA's  were  selected  according  to  two  criteria,  shelter/population  ratio 
and  population  density,  by  a  procedure  known  as  judgment  sampling  [Reference  2], 
Data  from  judgment  samples  only  suggest  or  indicate  conclusions,  which  is  precisely 
the  purpose  intended  in  this  analysis. 

Prior  to  selecting  the  SMSA's,  distribution  functions  of  population  density 
and  shelter/population  ratio  were  considered  (Figures  1  and  2).  The  distributions 
represent  200  of  the  213  SMSA's  listed  in  the  National  Location  Code  [Reference  3]. 
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Fig.  1.  Distribution  of  Population  Density  (200  SMSA's) 


A  .8  L2  L6  2,0  2.4  2.8  3.2  3.6  4.0  4.4  4.8  5.2  5.6 


Shelter/Population  Ratio  (spaces/person) 


|  Fig.  2.  Distribution  of  Shelter/Population  Ratio  (200  SMSA's) 


Two  constraints  were  placed  on  the  sampling  procedure: 

1.  SMSA's  with  population  exceeding  500,000  were  excluded.  This  constraint 
was  dictated  by  the  time  required  to  solve  the  algorithm  used  to  allocate 
people  to  shelter. 

2.  Multiple  SMSA's  (e. g. ,  Greensboro-High  Point,  N.  C.)  were  excluded.  In 
most  cases,  these  areas  require  special  consideration  because  of  complex 
movement  restrictions. 

Hie  shaded  portion  in  Figures  1  and  2  shows  the  distribution  (after  elimination 
of  multiple  SMSA's  and  SMSA's  with  population  greater  than  500,000)  frcm  which  the 
judgment  sample  was  taken.  The  characteristics  of  the  selected  SMSA's  are  listed  in 
Table  II,  Residential  Population  and  Fallout  Shelter  Characteristics  for  Eight  SMSA's. 

In  order  to  visualize  the  characteristics  of  the  selected  SMSA's  and  their 
relationships,  consider  the  joint  distribution  of  population  density  and  shelter/ 
population  ratio  as  shown  in  Figure  3,  Joint  Distribution  of  Population  Density  and 
Shelter /Population  Ratio  for  200  SMSA's.  If  we  consider  the  sample  to  represent  all 
SMSA's  lying  within  the  boundaries  formed  by  the  line  segments  connecting  each 
selected  SMSA,  then  146  of  the  200  will  be  represented.  (It  should  be  pointed  out, 
however,  that  all  SMSA's  within  the  boundaries  described  do  not  satisfy  the  popula¬ 
tion  constraint  of  being  less  than  500,000;  also  some  multiple  SMSA’s  are  included.) 

B.  Selection  of  Fallout  Environments 

One  form  of  available  data  providing  reference  intensity  (Iq)  and  associated 
arrival  times  for  fallout  (T^)  for  all  SMSA's  is  the  RISK  war  game  analysis 
[Reference  4],  In  order  to  use  representative  numbers  for  the  fallout  environment, 
data  from  an  0CD  study  in  a  form  similar  to  that  described  in  Reference  [4]  were  used 
to  derive  distribution  functions  for  both  I  and  T.  (Appendix  A).  Figures  4,  Cumu,- 
lative  Distribution  Function  for  Reference  Intensity,  and  5,  Cumulative  Distribution 
Function  for  Time  of  Arrival  of  Fallout,  relate  the  probability  of  occurrence  at  each_„. 
level  of  Iq  and  T^,  respectively,  for  the  0CD  study  data.  The  sample  size  (N)  was 
20,000;  i.e.,  100  attacks  on  each  of  200  SMSA's.  The  available  data  were  aggregated 
in  such  a  way  that  it  would  be  inaccurate  to  associate  a  joint  probability  of  occur¬ 
rence  of  Iq  and  T^.  However,  we  have  assumed  a  lower  limit  on  the  probability  of 
occurrence  of  a  given  combination  of  Iq  and  T^  to  be  the  product  of  their  individual 
probabilities.  By  assuming  complete  dependence,  we  can  arrive  at  some  upper  limit 
on  the  probability  of  occurrence  for  the  same  combination.  Thus,  we  can  describe 
some  probability  of  ~  terval  within  which  the  joint  probability  of  IQ  and  T^  must  lie. 


People 
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Joint  Distribution  of  Population  Density  and  Shelter/Population  Ratio  for  200  SMSA's 


niiuiiiim 


Cumulative  Distribution  Function  for  Reference  Intensity 


Time  of  Arrival,  (hours) 

Fig.  5.  Cumulative  Distribution  Function  for  Time  of  Arrival  of  Fallout 


For  example,  the  fallout  environment  will  be  I  <  3000  r/hr  and  T,  >  1  hour  between 

O  “  A  “ 

81  percent  and  90  percent  of  all  occurrences  (see  Table  III  for  method  of  computing). 

Six  fallout  environments  (combinations  of  and  Iq)  were  selected  and  used  to 
test  the  effectiveness  of  various  shelter  systems  (see  Table  III). 

TABLE  III 

Fallout  Environments  Msed  in  8  SMSA  Analyses 


I  T 

o  A  Approximate  Range  of 

(r/hr)  (hr)  Cumulative  Probability 


<  600 

£5.5 

367* 

- 

** 

607. 

2000 

- . .  • 

1.5 

- . . . - . 

727. 

857. 

3000 

1 

817. 

- 

907. 

*** 

5000 

1 

867. 

- 

957. 

*** 

10000 

1 

897. 

9S7. 

30000 

IHrk 

1 

.907. 

- 

99.97. 

Lower  Bound ;  Since  it  is  generally  true  that  late  arrival  times  correspond  to 
low  values  of  Ic»it  approximately  can  be  said  that: 

P  (ta  >  t  |  io  <  i)  >  P  (ta  >  t) 


so  that 

P  [Io  <  i,  TA  >  t)  -  P  (Io  <  i)  •  P  {Ta  £  t  |  IQ  <  i)  >  P  (Io  <  i)  •  P  (Ta  >  t) 
for  example: 

P  {I  <  600,  T.  £  5.5)  >  P  (I  <  600)  •  P  {T.  >  7}  £  (.60)  (.60)  £  .36 

Upper  Bound :  Since  a  probability  is  always  <1,  it  is  true  that: 

P  [Io  <  i,  TA  >  t)  -  P  (Io  <  i)  •  P  (Ta  >  t  [  ID  <  i)  <  P  (I0  <Tf  *1  ; 
and  similarly 

P  (Io  <  i,  ta  >  t)  <  1  •  P  (Ta  >  t) 
so  that 

P  {I0  <  600,  Ta  >5.5  }<  min  [p  {Iq  <  600 )  ,  P  (TA  >  7  ]  ]  <  .60 

In  some  cases  P  (TA  >  t)  was  not  known.  In  this  case,  the  uppe.  bound 
P  (IQ  <  i,  Ta  >  t)  £  P  (Iq  <  i)  should  be  used. 


-  10 


IV.  SM$A  ANALYSIS 


A. 


Casualty  Computations  and  Shelter  Utilization  Assumptions 

The  model  of  SMSA's  described  in  the  previous  sections  was  used  to  calcul 
the  number  of  casualties  in  each  SMSA  at  each  of  the  fallout  environments  desc 
in  Table  III.  (The  term  "casualty"  in  this  report  refers  to  anyone  receiving 
Maximum  ERD  greater  than  200r.)  Casualties  were  computed  for  each  of  the  four 
ment  and  shelter  allocation  rchemes  and  were  used  to  rank  their  relative  effec 
ness. 


ate 

ribed 

a 

move- 

Itive- 


Only  one  of  the  four  movement-to-shel ter  procedures  allowed  for  movement 
fallout.  In  each  of  the  four  procedures,  those  persons"  "not sheltered  in  NFSS 
(including  PF  Category  1)  were  gi^en  a  PF  of  2.  ^Residential  basements  were  no 
employed  in  these  analyses. 


Ln 

spaces 


It  is  recognized  that  residential  basements  afford  adequate  protection  in 
intensity  fallout  environments.  However,  excluding  these  spaces  will  only  affect 
the  casualty  computations  in  one  fallout  environment  of  our  analysis;  i.e.,  at 
reference  intensity  equal  to  600r/hr  and  time  of  arrival  equal  to  5.5  hours,  rtsiden 
tial  basements  afford  adequate  fallout  protection.  Three  SMSA's  of  the  sample 
shelter  in  residential  basements  for  more  than  107.  of  the  population — Pittsfield, 
Waterloo,  and  Richmond  [Reference  5].  Of  the  thre|,  only  Pittsfield  has  a  shelter 
deficit.  Therefore,  in  the  analysis  presented  here,  only  in  Pittsfield  could  the 
total  number  of  casualties  be  reduced  significantly  by  using  residential  basement 


low 


have 


data,  and  then  only  in  the  least  severe  fallout  environment. 


B. 


Analysis  of  the  8  SMSA's  1 

The  following  paragraphs  present  a  casualty  analysis  for  each  of  the  8  SM£>A's. 


They  are  presented  in  ascending  order  of  shelter  spaces  to  population  ratio.  Casu¬ 
alties  for  the  four  movement  assumptions  and  eight  fallout  environments  are  presented 
with  and  without  spaces  obtainable  in  NFSS  shelter  through  added  ventilation  or 
shielding.  An  analysis  is  made  for  each  SMSA  of  the  effect  on  casualties  of  the 
movement  assumptions. 


1.  Orlando.  Florida 


e  four 


The  SMSA  has  a  residential  population  of  263,540  and  a  population  density 

I 

of  4170  people  per  square  mile.  The  population  density  is  slightly  less  than 
the  average  of  the  200  SMSA's. 


I 
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Fig.  6.  Cumulative  Distribution  of  Population  and  Shelter 
in  Orlando,  Florida 

Approximately  20  percent  of  the  population  can  be  sheltered  in  spaces 
located  in  Phase  1  of  the  NFSS.  Only  5  percent  of  the  population  can  be 
sheltered  in  spaces  having  a  protection  factor  of  at  least  100  and  minimum 
volume  requirements  (500  cubic  feet  per  person  with  no  ventilation  or  10  square 
feet  with  adequate  ventilation).  _ _ 

The  plot  of  the  cumulative  distribution  of  population  and  shelter  as  a 
function  of  radial  distance  from  the  city  center  shown  in  Figure  6  indicates 
a  shelter  deficit  throughout  the  SMSA;  further,  almost  no  shelter  spaces  exist 
beyond  a  two-mile  radius  of  the  population  center.  Due  to  the  small  number  of 
available  shelter  spaces,  both  the  allocation  restricting  movement  to  two 
miles  and  the  constant  movement  speed  restriction  result  in  the  same  number  of 
casualties  as  that  if  movement  were  unrestricted. 

Casuilty  computations  are  significantly  affected  by  the  assumption  of 
movement  restricted  to  within  the  SL;  at  least  this  is  true  for  the  less 
severe  fallout  environments.  The  expected  number  of  casualties  at  each  fallout 
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environment  is  shown  in  Table  IV.  Figure  8,  Percentage  of  Casualties  for  Each 
Fallout  Environment  in  Orlando,  Florida,  shows  expected  casualties  from  Table 
IV  as  a  percentage  of  the  total  SMSA's  population. 


TABLE  IV 

Expected  Casualties  for  Each  Fallout 
Environment  in  Orlando.  Florida 


Fallout 

Environment 

Expected  Casualties  (7.) 

Reference 

Arrival 

Movement 

Movement 

Constant 

Intensity 

Time 

Restricted 

Restricted  — 

—  Movement  —  • 

-  -  -  Movement . . 

(r/hr) 

Oir) 

to  own  SL 

to  2  Mi. 

Speed 

Unrestricted 

600 

5.5 

94 

78 

78 

78 

2000 

1.5 

94 

78 

.  78 

78 

3000 

1 

96 

89 

89 

89 

5000 

1 

98 

93 

93 

93 

10000 

1 

98 

.98 

98 

98 

30000 

1 

99 

99 

99 

99 

Data  from  Phase  II  of  the  NFSS  survey  [Reference  6 ]  points  up  that 
approximately  7  percent  of  the  population  could  be  sheltered  ii  improvable 
spaces.  The  cost  of  improvement  would  be  $91,000,  or  approximately  $6.50  per 
space.  A  graphical  representation  of  the  shelter  status  is  presented  in 
Figure  7,  Percentage  of  Population  Sheltered  for  Each  PF  Category  in  Orlando, 
Florida. 

_ _ _It  is  quite  obvious  that  this  area  would  suffer  severely -in  any  fallout - 

environment.  Also,  only  slight  advantage  would  be  gained  by  improving  existing 
spaces.  The  only  solution  appears  to  be  an  accelerated  shelter  development 
program. 

2.  Tucson.  Arizona 

Tucson  has  a  residential  population  of  265,660,  about  the  same  as 
Orlando,  Florida.  The  population  density  of  3000  people  per  square  mile  is 
below  the  200  SMSA  average.  A  look  at  the  relative  location  of  population  and 
shelter  indicates  much  the  same  posture  as  Orlando  (Figure  9,  Cumulative  Distri¬ 
bution  of  Population  and  Shelter  in  Tucson,  Arizona). 
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Percentage  of  Casualties  Percentage  of  Population  Sheltered 


Fallout  Environment 


1  600  5.5 

2  2000  1. 5 

3  3000  1 

4  5000  1 

5  10000  1 

6  30000  1 


Casualties  Avoided  by 
Assuming  Improvable  Spaces 
are  Improved  and  Utilized. 


Excluding  the  SL  Movement 
Restriction 


Fig.  8.  Percentage  of  Casualties  for  Each  Fallout  Environment  in  Orlando. 
Florida  (Total  Population  -  263,540) 
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300 


Radial  Distance  From  Center  of  City  (miles) 


Fig.  9.  Cumulative  Distribution  of  Population  and  Shelter 
is  Tucson,  Arizona 

According  to  NFSS  Phase  1  data,  24  percent  of  the  population  can  be  shel¬ 
tered  in  spaces  with  protection  factors  ranging  from  20  to  1000.  Only  a  small 
number,  6  percent  of  the  population,  can  be  sheltered  in  spaces  with 
PF  )>  100  and  the  required  area  and  volume.  However,  improving  the  existing 
spaces  would  provide  shelter  for  an  additional  16  percent  with  protection 
factor  of  at  least  100.  The  cost  of  the  improvements  would  be  approximately 
$522,000,  or  $7.50  per  space.. 

Figure  10,  Percentage  of  Population  Sheltered  for  Each  PF  Category  in 
Tucson,  Arizona,  shows  the  shelter  status  of  existing  and  improvable  spaces. 
Although  the  fallout  shelter  status  in  Tucson  is  somewhat  better  than  Orlando, 
the  effect  of  movement  alternatives  is  the  same.  (See  Table  V,  Expected 
Casualties  for  Each  Fallout  Environment  in  Tucson.  Arizona.)  There  just  aren't 
enough  shelter  spaces  available,  even  if  improvement  alternatives  are  incorpo¬ 
rated  and  spaces  in  PF  Categories  below  4  are  accepted. 

Figure  11,  Percentage  of  Casualties  for  Each  Fallout  Environment  in 
Tucson,  Arizona,  again  illustrates  the  reduction  in  casualties  expected  from 
assuming  improvable  spaces  are  improved  and  utilized. 
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Percentage  of  Population 
Sheltered  in  Improvable  Spaces 


|~j  Percentage  of  Population 

Sheltered  in  Existing  Spaces 


PF  Category 


Percentage  of  Population  Sheltered  for  Each  PF  Category 
in  Tucson,  Arizona  (Total  Population  =  265,660) 


Percentage  of  Casualties 
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Fallout  Environment 
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Fig.  11.  lercentage  of  Casualties  for  Each  Fallout  Environment  in 
Tucson,  Arizona  (Total  Population  -  265,660) 
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TABLE  V 


I 

i 


Expected  Casualties  for  Each  Fallout 
Environment  in  Tucson,  Arizona 


Fallout  1 

Environment 

Expected 

Casualties  (%) 

Reference 

Arrival 

Movement 

Movement 

Constant 

Intensity 

Time 

Restricted 

Restricted 

Movement 

Movement 

(r/br) 

(hr) 

fo  own  SL 

to  2  Miles 

Speed 

Unrestricted 

600 

5.5 

89 

74 

74 

74 

2000 

1.5 

89 

74 

74 

74 

3000 

1 

94 

90 

90 

| 

90 

5000 

1 

97 

95 

95 

.. 95 . 

10000 

1 

98 

97 

97 

97 

30000 

1 

99 

98  * 

98 

98 

3.  Charleston,  South  Carolina 

This  SMSA  is  somewhat  atypical  in  that  the  population  density — 12,930 
people  per  square  mile--is  very  large  and  the  shelter  population  ratio  of  0.33 
is  rather  low.  The  general  trend  is  for  high  population  densities  to  be 
associated  with  a  relatively  high  shelter/population  ratio  (see  individual 
points  in  Figure  3).  Despite  the  combination  of  high  population  density  and 
low  shelter/population  ratio,  the  relative  location — with  respect  to  the  popu¬ 
lation  center — of  the  population  and  shelter  is  essentially  the  same  as  in 
Orlando  and  Tucson.  The  shelter  deficit  exists  throughout j the  SMSA  (Figure  12, 
Cumulative  Distribution  of  Population  and  Shelter  in  Charleston,  South  Carolina) 
As  might  be  expected,  except  for  the  SL  restriction,  there  is  no  variation  in 
the  number  of  expected  casualties  regardless  of  the  method  used  to  allocate 
people  to  shelters.  The  expected  casualties  are  shown  in  Table  VI,  Expected 
Casualties  for  Each  Fallout  Environment  in  Charleston,  South  Carolina  and 
Figure  14,  Percentage  of  Casualties  for  Each  Fallout  Environment  in  Charleston, 
South  Carolina. 

The  NFSS  indicates  that  approximately  33  percent  of  the  population  can  be 
sheltered  in  existing  spaces  (including  Category  1).  However,  as  in  Orlando 
and  Tucson,  the  great  majority  of  spaces  provide  a  protection  factor  of  less 
than  100.  Only  7  percent  of  the  population  in  Charleston  can  be  sheltered 
in  spaces  of  100  PF  or  better.  Adding  possible  improvements  would  provide 
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Cumulative  No.  (Thousands) 


300 


Radial  Distance  From  Center  of  City  (miles) 


Fig.  12.  Cumulative  Distribution  of  Population  and  Shelter  in 
Charleston,  South  Carolina  . . . 
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Fig.  13.  Percentage  of  Population  Sheltered  for  Each  PF  Category 

in  Charleston,  South  Carolina  (Total  Population  «  216,382) 
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Percentage  of  Casualties 


Casualties  Avoided  by 
Assuming  Improvable  Spaces 
are  Improved  and  Utilized 


Fallout  Environment 

No. 

itrjrnm 

TA_£hii 
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600 

5.5 
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2000 
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1 
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30000 

1 

* 

it 

Excluding  the  SL  Movement 
Restriction 


Fig.  14.  Percentage  of  Casualties  for  Each  Fallout  Environment 

in  Charleston,  South  Carolina  (Total  Population  *  216,382) 
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TABLE  VI 


Expected  Casualties  for  Each  Fallout 
Environment  in  Charleston,  South  Carolina 


Fallout  Environment 

Expected  Casualties  (7.) 

Reference 
Intensity 
.  (r/hr) 

Arrival 

Time 

(hr) 

Movement 
Restricted 
to  own  SL 

Movement 
Restricted 
to  2  Miles 

Constant 

Movement 

Speed 

Movement 

Unrestricted 

600 

5.5 

83 

70 

70 

70 

2000 

1.5 

83 

70 

70 

70 

3000 

1 

92 

89 

89 

89 

5000 

1 

97 

95 

95 

95 

10000 

1 

99 

98 

98 

98 

30000 

1 

99 

99 

99 

99 

shelter 

for  another 

6  percent  of 

the ' population. 

The  total  cost 

of  improve- 

ments  would  be  $130,000  or  $8.44  per  space. 

4.  Mobile.  Alabama 

Mobile  has  a  residential  population  of  273,942,  which  is  the  second  largest 
of  the  eight  SMSA's  being  considered.  The  population  density  is  the  lowest  at 
1320  people  per  square  mile.  Although  the  SMSA  has  an  overall  shelter  deficit 
of  54  percent,  it  is  different  from  the  previous  three  SMSA's  in  that  a  shelter 
surplus  does  exist  within  a  two-mile  radius  of  the  population  center  (Figure  15, 
Cumulative  Distribution  of  Population  and  Shelter  in  Mobile,  Alabama).  In 
fact,  the  area  outside  the  two-mile  radius  has  shelter  for  only  8  percent  of 
the  population,  and  62  percent  of  the  total  population  is  located  in  the  same 
region.  Therefore,  about  6  percent  of  the  population  must  '  -s  moved  "downtown" 
to  the  available  shelter  spaces.  The  expected  numbers  of  casualties  for  the 
various  movement  restrictions  are  shown  in  Table  VII,  Expected  Casualties  for 
Each  Fallout  Environment  in  Mobile.  Alabama.  As  can  be  seen,  there  is  little 
variation  in  casualties  according  to  the  method  of  allocating  people  to  shelter, 
with  the  exception  of  the  extreme  case  where  the  population  is  restricted  to  its 
own  standard  location. 

Using  spaces  located  by  the  NFSS,  48  percent  of  the  population  can  be 
sheltered  in  spaces  with  PF  ranging  from  20  to  1000.  (Figure  16,  Percentage 
of  Population  Sheltered  for  Each  PF  Category  in  Mobile,  Alabama).  Approximately 


0 


Total  Population  ■  273,942 


Percentage  of  Population  Sheltered 


Fig.  16.  Percentage  of  Population  for  Each  PF  Category  in  Mobile 
Alabama  (Total  Population  -  273,942)  * 


-  24 


Percentage  of  Casualties 


TABLE  VII 


Expected  Casualties  for  Each  Fallout 
Environment  in  Mobile,  Alabama 


Fallout 

environment 

Expected  Casualties  (7.) 

Reference 

Arrival 

Movement 

Movement 

Constant 

Intensity 

Time 

Restricted 

Restricted 

Movement 

Mbvemen  t 

(r/hr) 

(hr) 

to  own  SL 

to  2  Miles 

Speed 

Unrestricted 

600 

5.5 

78 

. . _  .  54  _ 

— .  54 

54 

2000 

1.5 

78 

54 

54 

54 

3000 

1 

91 

83 

83 

83 

5000 

1 

96 

.  93 

93 

93 

10000 

1 

97 

96 

96 

_ _  96  _ 

30000 

i 

99 

99 

99 

99 

15  percent  of  the  population  can  be  sheltered  in  spaces  having  FF  of  at  least 
100.  Improving  existing  spaces  would  provide  shelter  for  an  additional  9  per¬ 
cent  of  the  population  at  a  cost  of  $170,000  or  $6.90  per  space. 

5.  Macon.  Georgia 

Macon  has  very  nearly  the  same  shelter /population  ratio,  0.50,  as  Mobile; 
however,  the  population  density  is  much  higher.  The  total  population  is  141,249, 
and  the  population  density  is  4650  people  per  square  mile.  Macon,  like  Mobile, 
has  a  shelter  surplus  near  the  population  center  of  the  SMSA  and  very  few 
spaces  outside  a  two-mile  radius  of  the  center  (Figure  18,  Cumulative  Distribu¬ 
tion  of  Population  and  Shelter  in  Macon,  Georgia). 

According  to  NFSS  data,  50  percent  of  the  population  can  be  sheltered,  and 
20  percent  of  the  population  would  have  at  least  a  PF  of  100.  An  additional 
14  percent  of  the  population  could  be  sheltered  in  improvable  spaces.  Ihe  cost 
of  the  improvements  would  be  $94,000  or  $4.75  per  space. 

The  distribution  of  shelter  spaces  by  PF  category  shows  an  increase  in  the 
fraction  of  spaces  having  a  PF  of  100  or  better  over  SMSA's  previously  described 
(Figure  19,  Percentage  of  Population  Sheltered  for  Each  PF  Category  in  Macon, 
Georgia).  As  was  the  case  for  all  previously  examined  SMSA's,  the  largest 
increase  in  casualties,  for  all  movement  restrictions — as  shown  in  Table  VIII, 
Expected  Casualties  for  Each  Fallout  Environment  in  Macon.  Georgla--occurs  when 
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Radial  Distance  From  Center  of  City  (miles) 

Fig.  18.  Cumulative  Distribution  of  Population  and  Shelter  in 
Macon,  Georgia 


Percentage  of  Population 
Sheltered  in  Improvable  Spaces 


PF  Category 

Fig.  19.  Percentage  of  Population  Sheltered  for  Each  PF  Category 
in  Macon,  Georgia  (Total  Population  -  141,249) 
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Fig.  20.  Percentage  of  Casualties  for  Each  Fallout  Environment 
in  Macon,  Georgia  (Total  Population  ■  141,249) 


TABLE  VIII 


EXP 


ected  Casualties  for  Each  Fallout 
Environment  In  Macon,  Georgia 


Fallout  Environment 

Expected  Casualties  (7°) 

Reference 

Arrival 

Movement 

Movement 

Constant 

Intensity 

Time 

Restricted 

Restricted 

Movement 

Movement 

(r/hr) 

(hr) 

to  own  SL 

to  2  Miles 

Speed 

Unrestricted 

600 

5.5 

77 

46 

46 

46 

2000 

1.5 

77 

46 

46 

46 

3000 

1 

85 

73 

73 

73 

5000 

1 

87 

84 

84 

84 

10000 

1 

90 

86 

86 

86 

30000 

1 

97 

97 

97 

97 

the  fallout  reference 

intensity  changes  from  2000 

r/hi  to  3000  r 

/hr.  The  reason 

for  the  sharp  increas 

e  is  due  to 

i 

the 

ineffectiveness  of  Category 

1  spaces  at 

3000  r/hr. 

Of  the  to 

cal  number  of  spaces  available,  more  than  half  of  them  are 

Category  1 

spaces.  A 

characteristic 

i 

of  Macon  not 

encountered  in 

the  previously 

considered 

SMSA's  Is 

the  increase 

in 

casualties  at 

the  30,000  r/hr  fallout 

reference 

intensity. 

The  SMSA's, 

previously  considered,  experienced  rather 

constant  marginal  increases  over 

the 

range  from  3000  r/hr  to  30, 

000  r/hr. 

Macon,  however,  has  a 

relatively  large  number  of  Category  5  and 

6  spaces  (enough 

for  12  percent  of  the 'population)  which  are  ineffective  at  reference  intensity 
of  30,000  r/hr,  thereby  causing  the  increase  in  casualties  when  Iq  »  30,000  r/hr 

6.  Pittsfield.  Massachusetts 

Pittsfield  has  the  smallest  population  in  the  sample  (73,839),  and  the 
second  smallest  population  density  (1420  people/square  mile);  however,  its 
shelter /population  ratio  is  the  third  largest  (0.66).  A  shelter  surplus  exists 
within  a  one-mile  radius  of  the  population  center  (Figure  21,  Cumulative  Distri¬ 
bution  of  Population  And  Shelter  in  Pittsfield,  Massachusetts).  It  is  also 
the  only  SMSA  with  enough  improvable  spaces  to  erase  the  shelter  deficit.  By 
using  the  improvable  spaces,  47  percent  of  the  population  could  be  sheltered 
in  spaces  having  a  PF  of  at  least  100.  (Figure  22,  Percentage  of  Population 
Sheltered  for  Each  PF  Category  in  Pittsfield,  Massachusetts), 


lative  No.  (Thousands) 


Radial  Distance  From  Center  of  City  (miles) 

Fig.  21.  Cumulative  Distribution  of  Population  and 
Shelter  in  Pittsfield,  Massachusetts 
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Fig.  23.  Percentage  of  Casualties*  for  Each  Fallout  Environment 

in  Pittsfield,  Massachusetts  (Total  Population  ■  73,839) 
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As  in  all  the  other  SMSA's,  intense  fallout  environments  produce  a  high 
casualty  rate  as  shown  in  Table  IX. 


TABLE  IX 

Expected  Casualties  for  Each  Fallout 
Environment  in  Pittsfield.  Mass. 


Fallout 

Environment 

Expected  Casualties  (%) 

Reference 

Arrival 

Movement 

Movement 

Constant 

Intensity 

Time 

Restricted 

Restricted 

Movement 

Movement 

(r/hr) 

(hr) 

to  own  SL 

Co  2  Miles 

Speed 

Unrestricted 

600 

5.5 

64 

43 

43 

43 

2000 

1.5 

64 

43 

43 

43 

3000 

1 

72 

66 

66 

66 

5000 

1 

81 

81 

81 

81 

10000 

1 

87 

87 

87 

87 

30000 

1 

.  90 

90 

90 

90 

The  variation  in  expected  casualties  is  quite  large  when  comparing  condi¬ 
tions  where  some  movement  is  allowed  and  those  conditions  where  movement  is 
restricted  to  the  SL.  For  example,  in  the  fallout  environment  which  would  be 
exceeded  in  only  about  15-30  percent  of  the  attacks  (i.e.,  Iq  ■  2000  r/hr, 

T^  «  1.5  hrs),  approximately  64  percent  of  the  population  become  casualties 
if  movement  is  restricted  to  within  the  SL  while  only  43  percent  are  casualties 
with  any  other  mo\ement  restriction.  The  point  being  made  is  that  casualty 
computations  are  sensitive  to  movement  restrictions  in  the  majority  of  expected 
"fallout  environments. 


Thus  far  the  sensitivity  has  been  evident  only  when  movement  is  restricted 
to  the  SL  compared  to  movement  outside  the  jSL.  In  the  remaining  SMSA's, 
casualties  becomt  sensitive  to  all  four  movement  restrictions. 


7.  Waterloo.  Iowa 

Waterloo  is  characterized  by  a  shelter  surplus  over  the  entire  SMSA 
(Figure  24,  Cumulative  Distribution  of  Population  t.nd  Shelter  in  Waterloo,  Iowa). 
The  characteristics  of  Waterloo  are:  population--122,482 ;  population  density-- 
2120;  shelter /population  ratio--1.10.  Being  fortunate  enough  to  have  enough 
spaces,  hovever,  also  creates  new  problems.  In  all  the  previously  discussed 
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Fallout  Environment 

D  No  Movement  Restriction,  No  Movement  in  Fallout. 

®  Casualties  Added— Movement  Restricted  to  2  MPH,  Movement  in  Fallout 
is  Allowed. 

■  Casualties  Added— Movement  Restricted  to  2  Miles,  No  Movement  in 
Fallout. 


S  Casualties  Added— Movement  Restricted  to  Own  SL,  No  Movement  in  Fallout. 


Fig.  25.  Sensitivity  of  Casualty  Calculations  to  Movement  Restr ictions— 
Waterloo,  Iowa  (Total  Population  -  122,482) 
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Fig.  27. 
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in  Waterloo, 


Casualties  for  Each  Fallout  Environment 
Iowa  (Total  Population  -  122,482) 


SMSA's,  Che  shelter  deficit  was  such  that  there  was  no  serious  problem  in 
utilising  all  the  spaces.  In  fact,  the  available  shelter  could  be  treated  as 
if  the  population  has  unrestricted  movement  within  the  SMSA.  In  Waterloo,  with 
its  shelter  surplus,  the  expected  number  of  casualties  is  quite  sensitive  to 
any  change  in  movement  restrictions.  For  example,  consider  the  results  of  the 
allocation  model  as  shown  in  Table  X.  Only  in  the  most  severe  fallout  environ¬ 
ments  are  the  casualty  calculations  insensitive  to  the  movement  restrictions. 


TABLE  X 


Expected  Casualties  for  Each  Fallout 
Environment  in  Waterloo.  Iowa 


Fallout 

environment 

Expected  Casualties  (7.) 

Reference 

Arrival 

Movement 

Movement 

Constant 

Intensity 

Time 

Restricted 

Restricted 

Movement 

Movement 

(r/hr) 

(hr) 

to  own  SL 

to  2  Miles 

Speed 

Unrestricted 

600 

5.5 

83 

49 

0 

0 

2000 

1.5 

83 

49 

25 

0 

3000 

1 

87 

54 

42 

22 

5000 

1 

89 

64 

58 

57 

10000 

1 

91 

68 

68 

68 

30000 

1 

96 

83 

83 

83 

These  same  results  are  shown  graphically  in  Figure  25,  Sensitivity  of 
Casualty  Calculations  to  Movement  Restrictions--V'aterloo,  Iowa.  Quite  obviously, 
some  movement  is  necessary  in  order  to  utilize  shelters  and  avoid  large  numbers 
of  casualties.  As  will  be  shown  later  (in  Section  V  and  Figure  31),  because 
of  the  gradual  buildup  of  fallout  particles,  some  movement  in  a  fallout  field 

is  permissible  for  reference  intensities  as  high  as  5000  r /hr, — Restricting - 

movement  seriously  degrades  the  effectiveness  of  existing  shelter  spaces. 

8.  Richmond.  Virginia 

This  SMSA  has  the  largest  population  (408,494)  of  all  of  the  SMSA’s  in  the 
Sample.  It  also  has  a  60  percent  shelter  surplus,  (Figure  28,  Cumulative  Distri¬ 
bution  of  Population  and  Shelter  in  Richmond,  Virginia);  and  by  incorporating 
shelter  improvements ,  approximately  60  percent  of  the  population  could  be  shel¬ 
tered  in  spaces  having  PF  of  at  least  100.  Here,  as  in  Waterloo,  casualty 
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Cumulative  No.  (Thousands) 


Radial  Distance  From  Center  of  City  (miles) 


Fig.  28.  Cumulative  Distribution  of  Population  and  Shelter  in 
Richmond,  Virginia 
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Fig.  30.  Percentage  of  Casualties  for  Each  Fallout  Environment 
in  Richmond,  Virginia  (Total  Population  *  408,494) 
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calculations  are  sensitive  to  proper  shelter  utilisation.  For  example,  consider 

a  fallout  environment,  I  ■  3000  r/hr,  T.  *  1  hour.  If  movement  is.  restricted 

o  A 

to  within  standard  locations,  the  expected  casualties  will  be  81  percent.  How¬ 
ever,  if  movement  is  unrestricted  (i.e.,  assume  strategic  warning)  there  will 
be  one  percent  casualties.  Obviously,  these  are  extreme  examples.  Using  an 
optimum  allocation  plan,  such  that  casualties  are  minimized  subject  to  a  move¬ 
ment  speed  restriction,  would  result  in  approximately  38  percent  casualties. 
Using  the  two-mile  restriction,  which  is  similar  to  the  optimal  plan  except 
that  it  allows  no  movement  in  fallout,  would  result  in  68  percent  casualties-- 
a  30  percent  increase.  The  significance  of  limited  movement  in  fallout  is  quite 
evident.  The  expected  casualties  at  various  other  fallout  environments  are 
shown  in  Table  XI. 


TABLE  XI 


Expected  Casualties  for  Each  Fallout 
Environment  in  Richmond,  Virginia 


Fallout  Environment 

Expected 

Casualties  (%) 

Reference 

Arrival 

Movement 

Movement 

Constant 

Intensity 

Time 

Restricted 

Restricted 

Movement 

Movement 

(r/hr) 

(hr) 

to  own  SL 

to  2  Miles 

Speed 

Unrestricted 

600 

5.5 

66 

47 

0 

0 

2000 

1.5 

66 

47 

36 

0 

3000 

1 

81 

68 

38 

1 

5000 

1 

88 

81 

50 

45 

10000 

1 

90 

88 

87 

70 

30000  • 

1 

93 

91 

90 

90 

V.  DATA  ANALYSIS 


V.  The  Relationship  Between  Time  in  Fallout  and  Protection  Factors 

An  analysis  of  ERD  shows  interesting  trade-offs  between  time  in  fallout  and 
:he  protection  factor  of  the  shelter  that  is  entered.  Even  when  the  reference  in- 
;ensity  is  5000  r/hr  and  the  time  of  fallout  arrival  is  only  one  hour,  considerable 
:ime  in  fallout  can  be  compensated  for  by  a  relatively  small  increase  in  the  protec- 
:ion  factor  of  the  shelter  eventually  entered.  For  example,  a  person  traveling  15 
ulnutes  in  such  a  fallout  field  (at  2  mph)  and  reaching  a  fallout  shelter  with  a  PF 
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of  91  would  sustain  no  more  dose  than  a  person  who  was  in  a  shelter  with  a  PF  of 
70  at  the  time  fallout  commenced.  Of  course,  this  type  of  comparison  considers,  only 
gamma  radiation.  Other  types  of  radiation  might  also  have  to  be  taken  into  considera¬ 
tion,  but  their  effect  should  be  negligible  if  inhalation  of  particles  is  avoided. 

Figure  31,  Trade-Off  Between  PF  and  Time  in  Fallout,  shows  the  trade-offs  of 
PF  and  time  in  fallout  for  the  cases  where  the  reference  intensities  are  5000  r/hr, 
_3000  r/hr ^  and  2000  r/hr.  The  latter  has  an  associated  time  of  arrival  of  1.5  hours; 
the  others  have  a  time  of  arrival  of  one  hour.  The  curves  are  based  on  a  buildup 
time  of  twice  the  time  of  arrival  and  a  linear  rate  of  intensity  buildup. 

In  the  case  of  reference  intensity  of  10,000  r/hr,  there  is  little  opportunity 
for  trade-off  between  PF  and  travel  in  fallout.  A  person 

environment  would  receive  a  dose  greater  than  200r  ERD  if 

% 

as  ten  minutes,  regardless  of  the  type  of  shelter  he  even 


unsheltered  in  such  an 
h^e  were  outside  as  much 
ually  entered. 


Relationship  Between  Reference  Intensity,  Protection  Factor,  and  Maximum  ERD 


The  ERD  calculations  also  provide  a  means  for  comparing  reference  intensity, 
protection  factor,  and  Maximum  ERD.  Figure  32,  Trade-Offs  Between  PF  and  Fallout 
Environment,  shows  these  relationships  for  the  reference  intensities  considered 
in  this  volume.  The  dashed  lines  show  the  effect  of  times  of  arrival  of  more  than 
one  hour. 

'  » 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  Discussion  j 

Shelter  utilization  analyses  necessarily  include  assumptions  concerning 
variables  that  are  difficult  to  measure.  These  variables  affect  utilization.  Among 
the  more  significant  variables  are:  . . . 

1.  Knowledge  on  the  part  of  people  about  the  shelter  which  they  should  use 
(shelter  assignment  plan,  discipline,  crowd  control,  etc.). 

2.  Time  required  to  reach  shelter  (warning  time,  start-up  time,  travel  time, 
searching  time,  interference  of  traffic  flow,  etc.). 

3.  Distance  to  shelter  (influences  both  time  and  knowledge,  for  the  knowledge 
may  be  less  if  distances  are  long). 

In  a  complete  analysis,  relations  among  these  three  and  other  variables  must 
be  derived  (for  they  are  obviously  not  independent),  and  these  in  turn  related  to 
the  shelter  utilization  pattern. 
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In  the  analyses  reported  here,  knowledge,  time,  and  distance  to  the  shelter 
were  introduced  in  a  semi-quanti^ative  way  by  varying  "movement  restrictions"  as 
follows: 

1 .  Movement  Restricted  to  Own  SL 

Knowledge  very  great,  and  time  and  distance  allowance  great  enough  to 
allow  best  shelter  utilization  within  one's  own  SL. 

2 .  Movement  Restricted  to  2  Miles  From  Own  SL 

Knowledge  very  great,  and  time  and  distance  allowance  great  enough  to 
allow  best  utilization  of  any  shelter  within  one's  Standard  Location  (SL) ,  or 
a  SL  not  more  than  two  miles  away  (rectangular  distance). 

3.  Movement  Unrestricted  Within  SMSA 

Knowledge  very  great,  and  time  allowance  great  enough  Lu  allow  best  shel¬ 
ter  utilization,  at  long  distances  anywhere  in  the  SMSA. 

4.  Optimum  Shelter  Utilization  (Transportation  Algorithm) 

Knowledge  exact,  to  allow  optimal  shelter  utilization  (minimum  casualties) 
anywhere  within  the  SMSA,  including  a  prior  calculation  of  dose  received  during 
transit  in  fallout  and  in  the  primary  shelter,  all  before  cr  during  fallout 
arrival.  Speed  was  restricted  to  2  tnph. 

We  have  a  mathematically  optimum  solution  in  Movement  Restriction  4.  A  similar 
ihelter  utilization  pattern  could  be  generated  by  other  assumptions,  operationally 
tore  realistic--as  seen  In  Table  XII,  Summary  of  Casualty  Estimated  (%  of  Population) 
'or  Reference  Intensity  of  2000  r/hr  and  Time  of  Arrival  of  1,5  Hours--for  unrestricted 
lovement,  for  example.  For  more  severe  environments,  the  two  casualty  figures  (hence 
helter  utilization)  become  identical.  This  case  and  unrestricted  movement  offer  a 
■enchmark  of  ideality  against  which  to  measure  other  shelter  utilization  patterns. 

.  Major  Conclusions 

Efficient  shelter  utilization  appears  to  be  a  critical  problem  when  the  shelter 
paces  exceed  population;  inefficient  utilization  can  lead  to  a  large  number  of 
voidable  casualties.  This  was  concluded  from  the  observation  of  the  rapid  increase 
n  casualties  for  a  specified  attack  environment  as  one  restricts  movement  progres- 
ively  from  the  SMSA  down  to  the  SL  in  which  the  shelteree  finds  himself  at  the  time 
f  warning.  For  the  most  shelter-abundant  cities  (Waterloo  and  Richmond  in  Table 
II),  the  percentage  of  casualties  under  a  2000  r/hr  fallouc  environment  is 
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progressively  reduced  from  a  tvo-city  average  ranging  between  757.  and  307.  for  limited 
movement  (SL  and  Constant  Movement  Restriction  in  Table  XII)  to  zero  for  movement 
unrestricted  within  the  SMSA. 


TABLE  XII 

Summary  of  Casualty  Estimates  (%  of  Population)  for 
Reference  Intensity  of  2000  r/hr  and  Time  of  Arrival  of  1.5  Hours 


SMSA 

Shelter/ 
Population 
Ratio  (%); 
Categories 

1  and 
Better 

Movement 

Restricted 

To  Own  SL; 

Movement 

Completed 

Before 

Fallout 

Constant  Movement 
Speed  of  2  mph 
Restriction 
Movement  Permitted 
in  Fallout  (Trans¬ 
portation 
Algorithm) 

Movement 
Unrestricted 
Within  SMSA; 
Movement 
Completed 
Before 
Fallout 

Orlando 

20 

94 

78 

78 

Tucson 

24 

89 

74 

74 

Charleston 

33 

83 

70 

70 

Mobile 

48 

78 

54 

54 

Macon 

50 

77 

46 

46 

Pittsfield 

66 

64 

43 

43 

Waterloo 

110 

83 

25 

0 

Richmond 

160 

66 

36 

0 

For  the  most  shelter-poor  cities  (Orlando  and  Tucson),  the  corresponding  reduc¬ 
tion  is  from  a  two-city  average  of  917.  casualties  to  a  two-city  average  of  767.  casu¬ 
alties,  a  feasible  savings  of  15%  of  the  population  who  would  otherwise  be  casualties 
This  is  to  be  compared  to  the  757.  savings  to  the  more  shelter-abundant  cities.  The 
lack  of  sensitivity  to  movement  restrictions  in  shelter-poor  cities  generally  is  due 
Co  the  fact  that  there  are  enough  people  in  any  standard  location  to  completely 
utilize  available  shelter  in  that  SL.  Thus,  even  the  most  restrictive  shelter  utili¬ 
zation  constraint  leads  to  good  utilization. 

Within  broad  limits,  tne  higher  the  ratio  of  shelter  is  to  population,  the 
greater  is  the  payoff  in  casualties  saved  by  shelter  utilization  planning. 

Insofar  as  movement  distance  restrictions  are  derived  from  time  available  for 
sheltering,  one  can  deduce  the  substantial  payoff  potential  in  effective  warning 
means,  shelter  planning,  and  indoctrination  of  the  population--particularly  in 
shelter-abundant  areas. 

C.  Secondary  Conclusions 

1.  Movement  in  fallout  to  a  better  shelter  theoretically  is  warranted  in 

many  circumstances.  The  quantitative  trade-offs  between  movement  time  in 
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fallout  and  protection  factor  of  the  shelter  are  shown  in  Figure  31.  Un¬ 
feasibility  of  applying  this  principle  in  the  transattack  environment  is 
questionable,  however.  It  is  a  special  case  of  remedial  movement,  but  with  a 
very  short  decision  time  and  a  changing  environment. 

2.  Neither  city  size  nor  population  density  (one  of  the  city  selection'  criteria) 
which  might  be  supposed  to  have  an  effect  on  juxtaposition  of  people  arid  shelter, 
could  be  shown  to  be  significant  in  this  analysis.  This  is  concluded  from 
examining  the  difference  in  casualties  in  each  attack  environment,  between 
Movement  Restriction  4,  (movement  in  fallout,  but  limited  movement  speed),  and 
Restriction  3  (complete  freedom)  as  a  function  of  city  size  and  population 
density.  Small  effects  which  likely  exist  either  are  swamped  by  the  more 
critical  shel ter/population  ratio,  or  are  undetectable  from  the  "noise"  of 
variable  PF  distributions  in  cities  otherwise  similar  (size  and  population 
density).  The  casualties  differ  between  Restrictions  3  and  4  only  for  Waterloo 
(low  density,  small  size)  and  Richmond  (high  density,  large  size)  at  the 
intermediate  attack  environments. 

I).  Recommcndat  ions 

This  analysis  has  shown  the  critical  sensitivity  of  casualty  estimates  to 
shelter  utilization,  and  therefore,  to  the  corresponding  movementrto-shelter 
iissumptions ,  particularly  for  shel ter -abundant  areas.  Accordingly,  further  research 
is  recommended  to  give  a  more  quantitative  insight  into  the  problem  of  movement -to- 
shelter  (time,  training  and  discipline,  and  distance  to  move).  Better  understanding 
can  lead  to'  more  effective  operational  plans  and  programs  in  terms  of  preattack 
shelter  assignment  planning,  warning  requirements,  training,  etc. 
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Appendix  A 

Development  of  Probabilistic  Fallout  Environments 

1.  INTRODUCTION 

A.  Cumulative  Distributions  Mathematically  Defined 

Given  a  simulated  nuclear  attack,  the  pattern  of  the  attack  is  used  to  compute, 
on  a  probability  basis,  the  severity  of  the  fallout  environment  at  each  of  200  SMSA's 
throughout  the  U.  S.  One-hundred  attack  patterns  comprise  one  RISK  attack.  Thus  for 
one  RISK  attack,  each  SMSA  will  be  subjected  to  100  possible  fallout  environments. 

Let  I  and  T^  dnote  the  Reference  Intensity  and  Time  of  Arrival,  respectively. 
Each  component  is  divided  into  subgroups: 

A  similar  distribution  is  derived  for  time  of  arrival  T. : 


i 

I 

0 

ta 

(r/hr) 

(hr) 

1 

0-100 

0-1 

2 

101-300 

1-2 

3 

301-1000 

2-3 

A 

1001-3000 

3-4 

5 

3001-10000 

4-8 

6 

10001-30000 

8-12 

If  X.  .  denotes  the  number  pf  times  I  occurs  in  the  interval  for  the  i1*1 

ij  0  th 

SMSA  then  the  probability  of  occurrence  for  the  1  interval  is  estimated  by 
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where  f  represents  the  number  of  attack  patterns.  Similarly,  if  denotes  the 
number  of  occurrences  of  in  the  i  interval  for  the  SMSA,  we  have 


A- 1 


6 


P(Yt) 


2  "Vi 
>;  \ 

x± 


ij 


200 


f.  i-1. 


(A-2) 


The  discrete  cumulative  distribution  function  for  each  variable  is  developed  by 
summing  the  individual  probabilities. 
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Tin.'  Populat  ion  to  Shelter  Algorithm 


I 

I.  INTRODUCTION 

■I 

A.  Expression  of  Content  .  /' - - -  • . . . . 

This  appendix  describes  in  detail  the  model  used  to  estimate  casualties  within 
an  SMSA.  Also  included  is  the  computer  program  used  in  calculating  the  Maximum  ERD 
associated  with  movement  between  standard  locations. 

1.  Development  of  a  Mathematical  Model  to  Calculate  Maximum  Equivalent 
Residual  Dose  (Stage  1) 

The  Model  and  Program  described  in  this  section  is  a  modified  version  of 
the  Mainline  Model  and  Program  discussed  in  Appendix  of  Volume  III,  "A 

Generalized  Sensitivity  Analysis  of  CD  Systems."  The  greatest  difference  between 
this  algorithm  and  the  Mainline  Model  of  Appendix  B  of  Volume  III  is  in  the  input/ 
output  part  of  the  program.  By  necessity  each  being  oriented  to  a  different  re¬ 
search  objective. 

The  equation  describing  ERD  at  t  hours  is  the  sum  of  the  reparable  and 
irreparable  dose  over  finite  time  increments. 


Fig.  B-l .  Illustration  of  Finite  Different  Formulation 
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T.  *nd  T  denote  the  time  of  arrival  and  time  of  cessation  of  fallout 

Ac 

respectively. 

*  Hr 

Let  Dp  represent  the  reparable  dose,  and  the  irreparable  dose  for  the 
p**1  interval.  Then  . 


D*  -  [1  -  (a)  (b)  ]  D*  l  +  (1-f)  ^11^- 


(B-l) 


**  ** 

D„  -  Dn  +  (f)(r)(a) 
P  P"1 


(B-2) 


where 


(1)  a  is  the  constant  time  increment  in  hours. 

(2)  b  is  the  recovery  rate  (0.1%/hr.). 

(3)  f  is  the  nonreparable  part  of  the  dose. 

(4)  PF  is  the  protection  factor  for  the  shelter  used. 


r  (t  -  a/2  -  Ta) 


T.  <  t  <  T 
A  —  —  c 


(B-3) 


T  =  2T. 
c  A 


r  ■  a  t 


,  t  >  T 

'  f' 


(B-4) 


The  ERD  for  the  p  time  increment  is 


*  ** 

D  -  (D  +  D  )  I  , 
P  P  P  o 


(B-5) 


and  the  Maximum  ERD  occurs  when  D  >  D  . , 

P  P+1 


Ihe  previous  derivation  is  based  on  the  assumption  that  a  person  enters  a 
shelter  with  protection  PF  before  fallout  arrives;  i.e.,  if  the  time  required 
to  move  from  standard  location  i  to  j  is  T^ ^ ,  then  the  assumption  is  that 
T^j  <  Ta  .  If  Tf i  >  Ta  ,  then  in  Equation  B-l,  PF  must  equal  1  for  T{  ,  >  t  . 


lij  A 


When  T^  becomes  less  than  or  equal  to  t  ,  (T^  <  t)  the  appropriate  PF  is 
empLoyed  and  the  original  derivation  is  valid. 


The  computer  program  for  Stage  1  is  described  in  Sections  3,  4,  5  of  this 
Appendix. 
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Th«A«  s  i  z  — se-u  x,.x!«  \  (St  aze  2) 


The  assignment  n,>»e !  i «  t-a«e  1  on  a  linetr  ft  •, .raining  technique  called 
*. '  ,  h  i  •  »  M*  *  -  *  .  "  Itiis  technique  is  most  commonly  employed  in 

'.it  '  t*  *  l:.e  cost  of  shipping  a  homogeneous  product  manufactured  in  "n"  mill 
4*\ch  in  a  different  geographical  location,  and  shipped  to  consumers  at  "m" 
different  destinations.  The  analogy  is  made  that  people  represent  the  product 
located  in  "n'L standard  locations  and  fallout  shelter  spaces  represent  the 
warehouse  capacities  at  "m"  standard  locations.  The  cost  variable  is  the 
equivalent  residual  dose  received  if  one  moves  from  standard  location  i  to 
standard  location  j  where  j  has  a  protection  factor  of  k  .  ^ 

-  Symbolically,  the  model  formulation  is:  . "  -  * . I------  r 


Let  represent  the  ERD  defined  previously:  |'|J' 

sijk  ■  f<R-T*-v-sij-prk>  <b-m 

R  "»  fallout  reference  intensity  (r/hr) 

T.  —  arrival  time  of  fallout  (hr) 

A 

V  ~  movement  speed  (mi/hr) 

distance  from  SL  i  to  SL  j  (mi) 

PFk  ~  protection  factor. 

■  *' 

The  model  which  calculates  s,„  is  defined  in  Section  A.l.  Then  assume  the 

ijk 

ERD  obtained  in  moving  from  any  SL  i  (i  ■  1 ,  . . . ,  n)  to  any  SL  j  (j*>l,  .... 

✓ 

with  a  Drotection  factor  category  k  (k  ■  1,  8)  Is  known.  Thus,  we  have 

a  3-dimensional  matrix  (n  x  m  x  8)  with  each  element  representing  a. "cost" 

t  Hi  >'  . 

[or  Maximum  ERD.  To  simplify  the  matrix,  let  jk  represent  the  j  destination 
th ' 

and  k  protection  category.  What  this  means  is  that  each  protection  category 
is  treated  as  a  separate  and  distinct  destination,  where 


jk  -  8 ( J - 1 )  +  k  . 


(B-7) 


For  example,  let  the  cost  matrix  be  represented  by  D  and  let  s  represent 

the  Maximum  ERD  for  a  person  leaving  standard  location  3,  moving  to  standard 
location  5,  and  entering  a  protection  factor  category  8.  The  element  in  D 
corresponding  to  s3  5  g  is  where 


i  -  3 

jk  -  8(5-1)  +  8  -  40  . 


(B-8) 


There  fere 


*1  t  c  —  »  1  A  *nd  S  ,  -  D  . 

3,f  ,8  sr  3  .  -*0  3,-0 

To  summarize  briefly  the  derivation  of  the  cost  matrix: 

(a)  The  maximum  ERD  is  calculated  in  Stage  1  of  the  model  for  each 
element  of  a  3-dimensional  matrix  where  s  ^  represents  each  element 

— in  the  matrix.  ;;  ~  •'  . “""",7;  .  . . 

(b)  In  order  to  fit  the  2-dimensional  cost  matrix  scheme  of  the  transporta¬ 
tion  algorithm,  a  linear  transformation  was  performed  on  each  element 

whereby  s.  became  s,,,.v. 

ijk  Kjk) 


With  the  cost  matrix  in  the  desired  form,  the  constraint  equations  and 
objective  function  may  be  defined: 

(a)  The  population  allocated  from  origin  i  to  destination  jk(X  .  .) 

*  \  J  K/ 

must  not  exceed  the  initial  population  at  the  origin,  A^. 


Z  x  (jk)  <  A  ,  (i  «  1 . n) 

jk-1  1  1 


(B-9) 


Xi(jk)  -  0  * 

(b)  The  population  sheltered  at  destination  jk  must  not  exceed  the  shelter 
capacity  at  the  destination,  B.  ... 


Xi(  jk)  -  B(  jk)  *  Ok  “  1.  ....  8m) 


(B-10) 


Xi( jk)  -  0  * 


(c)  The  objective  function  to  be  minimized  is  the  number  of  casualties, 
subject  to  the  previously  defined  population  and  shelter  space  con¬ 


straints. 


n  8m 

min  z  «  Z  Z  c..,. .  x.  . . 

i-1  jk-1  l^k)  ^k) 


i  (  jk.) 


°’  di( jk)  -  200r 
l*  di(jk)  >  200r 


(B-ll) 
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ISL  Number  of  standard  !  >c  a  t  i  s  in  the  'jMSA 

LMAX  Maximum  number  of  PF  categories 

Z  Radiation  decay  constant 

R1  Reference  intensity 

TA  Arrival  tine  of  fallout 

VEL  Movement  speed 

B  Recovery  rate 

F  Recoverable  fraction 

DMAX  Maximum  acceptable  ERD 

PF1  Primary  PF 

PF2  Secondary  PF 

PF3  Equivalent  PF  in  third  shelter  period 

PF4  Equivalent  PF  in  fourth  shelter  period 

T1  Time  in  primary  shelter 

T2  Time  in  secondary  shelter 

T3  Time  in  third  shelter 

III  SMSA  code 

IP  Number  of  populated  standard  locations 

IS  Number  of  standard  locations  with  shelter  spaces 

IC  IS  x  IP 

RMAX  Peak  reference  intensity  of  buildup  function 

D1  Reparable  portion  of  total  dose 

D2  Irreparable  portion  of  total  dose 

DD  Equivalent  Residual  dose  (ERD) 

DP  Maximum  ERD 

TMM  Cumulative  time  since  time  of  arrival 

L . . . .  PF  category.  1  ~  ~  “ 
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I  BEGIN  LOOP 
TO  CALCULATE 
MAX.  ERD  FOR 
MOVEMENT  FROM 
3l(i)  TO  SL(j) 


CALC. 

<Jl)  DISTANCE 
TRAVELED 
(2)  TIME  IN 


Exit 
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Fortran  Listing 


76655  ISS9XI  1C56 

•  .  .  ..SC3NCV  SiM  .  .RH-QCD 

•  xe  r; 

SECTION  I  -  CALCULATION  OF  MO  VEMFn  T  TIMES 

OIMENSICN  XI 100 »  » VI  ICC) vn(l 01*01(10  I *02(101 ,CPIlC>,TMtlCI 
DIMENSION  PFiUO>,PF2(iOl  , S 1 10CC 1 , JR S AC  1 1 00 ) 

DIMENSION,  DDPIAG.AC.B) ,TFRm(AG,AC )  ,PflOO) 

rewinc  g 

I GNCRF  =C 

1  CONTINUE 
11=1 

R F AC  INPUT  TAPE  5 , 3C 1 , I $L , LM AX , ? 

IFUSl  171A.71A.2 

2  READ  INPUT  TAPE  5 , 5C 1 , R 1 , TA . VEL 
READ  INPUT  TAPE  5 , 5C2 ,B ,F ,DMAX 

1C  CONTINUE 

READ  INPUT  TAPE  5  ,  A  ,  PF 1 (II , PF 1 1  2  I , PF It  3  I , P F 1 ( A ) , PF 1 ( 5 » 

READ  INPUT  TAPE  5,A,PF1{6)  ,PF1C7),PF1I8)  ,PF1(9)«PF1UC» 

RE  AC  INPUT  TAPE  5»A»PF2I1),PF2C2)»PF2I3)»PF2(A),PF2(5I 
RF  AC  INPUT  TAPE  5  ,  A  ,  PF 2 1  6  I , PF2 C 7 ), PF 2  I  A ) ,PF2 1 9 ) ,PF2 1 1C » 
READ  INPUT  TAPE  5  ,  A ,  PF  3 ,  PF  A  ,  T 1 ,  T2,  T  3 
3C1  FORMAT I2I5*F1C. 2) 

5  C  1  FORMAT  I3F1C.2) 

5C2  FORMAT  I  3F1C.5) 

A  FORMAT  (5F1C.2) 

WHITE  CUTPUT  TAPE  6 , 8C  1  ,R 1 . T A , VFL , PF 1 f 1 ) 

8C1  FORMAT  (AF12.2) 

RE  AC  INPUT  TAPE  5,7C7 

707  FORMAT (ACHI  ) 

WRITE  OUTPUT  TAPE  9,7C7 
WRITE  CUTOUT  TAPE  6,707 
READ  INPUT  TAPE  5 , 7C5 . I II , I P , I S , 10 
WRITE  CUTPUT  T  APr  6 , 7C 6 , 1 P  ,  I S , IC 
WRITE  CUTPUT  TAPE  9 , 7C6 , 1 P  ,  I S , I C 
7C  5  FORMAT!  17, 1  13, 1 1C,  1 1G) 

?C6  FC  R_M  A_T  (  7HC  2  1 1 1 1  3  I  1  3  , 11 0 ,  I  1 C  1 
IF (  IGNORE ) 7061,706 1.7C62 
706  1  DC  703  M *  1 ,  I SL 

R E AC  INPUT  TAPE  5 , 70 1 , JR S AC ( M J , X { M ) , Y I M } 

701  FORMAT! 18, F12.C.F10.0)  ’ 

7C3  CONTINUE 
7062  CONTINUE 

RMAX  *  ( 2.0*T  A  !•* {-2 ) 

7C8_C0NTINUE  . 

WRITE  CUTPUT  TAPE  6,8CC 

flOC _ fCRMA_T  ( 5  5H1  I  J  PF  JJ  PRCB  MAX  ERD 

DO  305  1*1, ISL 

PC  .305  J*  I.,  ISO _  ..  .  .  _ 

5C9  DO  *  C 
Dll  *  C 
D22  *  C 
TMW  *  C  „ 

DPP  *  C 

_ t  *  ($cRTFnxm-x uj i”2 ) ^Rif m m - v i j i Lt-2ii/m_ .  ... 

T  *  T»0,062 
§90  CONTINUE 
5C A  CONTINUE 
5t£  .CONTINUE  ... 


TIME) 


1 


T  «  T-TA 

SECT  I CN  I  -  CALCULATION  OP  MOVEMENT  TIMES 

JF _ C T 1  5C_7»  5C8 1_508 

5C7  T  *  0. 

_ 5C8  CONTINUE _ 

IF(T)  599,599,510 

510  TMM  »  TMM»Q,  1 _ 

RM  *(RMAXMTMM-C.C5)/TA')*C.1 
IF { {RM»  10.  l-RMAX >515* 515,5 14 
514  l.-Z) 

__  _A2  a ( T  M  M-  .  1  }••{  l.-Z) 

RM«(A1-A2)/I l.-Z) 

- 515  01 1_»  y_l -.0-0. 1»B  1_* C 1 1  )♦  I  ( 1.0-F  > HI 

022  «  022* I F  *RM ) 

_ CC_  =  (011*022)»Rl 

IFIDD-CPP)  513,511,511  . 

_ 5ii  opp  =  cc 

IF! DPP- DM AX)513, 599,599 
513  I F  I  TMM- ( T- 0 .05)1  51C.599.599 
599  CONTINUE 

_C  _  SECT  I C  N_H  I 

5  CONTINUE 
_7  _  OC  101  L  =  l.LMAX 
OIL  I  *  CD 

_  DHL)  »  Dll  _ 

02  IL)  '=  02  2 
_DP(L)  *  OPP 
TM(L1  *  TMM 

_  IFICP(Ll-OMAX)  19,999,999 

15  I F ( TMM-T  A )  20,30,30 

_  TM{l)*C.5_ 

RML  =(RMAX*'(  ITMll  )-C.*25j/TA)  )«0.5‘ 

Ol(L)  *  (l.C-B/2.C)»01(U  +  (l.C-F)*RML/PFl(U 
02  I L  j  =  D2(L)*F»RML/PF1(L) 

OILI  =  R 1* I C 1 ( L ) *02  t  L  )  1 
IF  (OIL  )-DP  (D)  22,22,2 1 

21  TERM(I  ,J)  =  TMIL)  _  _ 

OP  (L)  *  DILI 

IF  I  CP  I  U-700.  122,999,999 
22  TAL  *  TA-0.25  I 

IF|TM(U-TAU2C,3C,30 
30  TMILI  =  TMIL)  +6.0 

_ M  *_T MIL  )  ••  I  1.0— Z1  _  _  _  < 

A2  s  ( f M I L )-6 . C ) *• I l.C-Z ) 

RML  *  (A1-A21/I l.C-Z) 

OKU  =  ( l.C-6.*B)*01  C  L  I  ♦  C  l.O-F)  •RML/PF1IL  ) 

.  02  I L  )  =  02(U*F*RML/PFl(U 

OIL)  *  R1MD1  (LKC2IL  )  ) 

_ LLl-SlU-DfJU  I,32i  32. 31. . 

31  TERM  I  I , J )  =  TM( L  ) 

_  OP(L)_  =  0(L  I  _ 

IF  (CPU  )-700. ) 32, 995,999 

_ U  IF(TM(L)-T1)3C,40,40 

40  CONTINUE 

4C1  TM(L)  »  TM  (L  )  ■♦■24.  C  _ 

A 1  *  TMIL)*»< l.C-Z) 

_ 1-6.Q)**  LUfl-i  )  .  . 

RML  -  (  A 1-A2 1/(1. C-Z ) 

I _ _ _ S££UflH_I  ~  CALCULATION  DF  MOVEMENT  TIN1S 


.  .  .  01 C L )  *  (l.£-2*.0»6)*01(M 

01(L)  =  01  (L  I  ♦  U  l.C-F  )*RML  >/PF?(U 
02(U  *  D2(L)*F*RML/PF2d) 

0(1)  =  Rl*(0l(L)*C2d)  ) 

IP  (OIL  )-OP(l)  )A2,A2,Al 
Al  TERM(I ,J)=TM(L) 

OP(L)  =  0 ( L  ) 

IF(CP(L)-700.  >A2,99< 9,991 
A 2  TAL  =  TUT2 

IF (TM(L J-TAl )  ACi.A9.A9 
A 9  IF(T3)999, 999,50 
5C  T*M  L  )  =  TM(L  )*168.0 
A 1  =  TM(t)»*(  l.C-Z) 

A2  =  (TMd)-6.0)*»(l.C-Z) 

RVL  =  (A1-A2)/(X.C-Z) 

Old)  =  (  1  .C-  16R.C-B  )*01  {  L  )♦  (  l.C-F  )  *RMl/PF3 
02(1)  =  02 (L ) ♦F*RML/PF3 
0  (  L  )  =  R1»(01(L)*C2(L)  ) 

IF(C(L  )-OP(U  )  52,52,51 
51  TERM ( 1 , J  j  =  TP { L ) 

OP(l)  =  D( L ) 

IF (CP(L)-700.) 52,999,999 
52  TAL  =  T1+T2+T3 

IF (TM(t)-TAL )5C»39,59 
5  9  CCnTINUF 

6Q  TM(l)  =  TM(U*720. 

A 1  *  TM(L)*«( l.C-Z) 

A2  *  (TM(L)-6.0)**(1.C-Z) 

Rml  =  ( A 1-A2 ) / ( l • C-Z ) 

Ol(L)  =  (  1. 0-720.  C»B)*C1(UM1.C-F)*R'*L/PFA 
1)2  (U  *  02  (L  )  ♦F»RML/PFA 
0(L)  =  Rl*(Cl(L)+D2d)  ) 

I F ( 0 ( L  )-OP(L )  ) 62 , 62 . 6 1 
61  TERM  ( I »J)=TM(t) 

DP(L)  =  OIL) 

I F ( DP ( L )-70C • 162,999,999 
6  2  CONTINUE 
GC  TG  6G 

999  OCP ( I , J,L )  =0P J  L ) 

1C1  CONTINUE 
3C5  CONTINUE 

OC  90  J* 1 , 1 SL 
_  DC.  90  J  =  Jjl5l 

or  90  l=T,Cmax 

OCP 1 1 , Jf L )  =  OOP ( J , I , L ) 

TFRM(I,J)  =  TFRM(J,I) 

9C  CONTINUE 

OC  91  1=1, I SL 

_ DC  -91  0=1,1  SL  .  .. 

DC  91  1=1, UMAX 
__  .  JJ=(J»LMAX)-dMAX-L) 

P~ROB=(CCP(  I  »J»L)-20C.  )/50C. 

_ I  F.1  PRC  £19  ’,9  A,  9  A 

93  NPRCB=C 

GC  TO  95 _ 

9 A  IF  (PRC 8-1.  )9A2»9A1,9A1 

SECTION  I  -  CALCULATION  OF  MOVEMENT  TIMES 


B-ll  - 


?AJ  NPRCB=1000 


GO  TO  <35 

942  NPRCB»PRQB»10CC. _ 

95  WRITE  OUTPUT  TAPE  6 , 709 ,  l‘,*J,  L «  J j'.NPROB  ,  OOP  Cl  ,~J  ,  L  5  ,  TERM  II  ,  J  ) 
WRITE  OUTPUT  TAPE  9 , 720 , I , J J , NPROR 


709 

FORMAT  1 1CH 

415, 110, 2F 1C. 3) 

72C 

FORMAT  J 10HB 

15,15,110) 

91 

CONTINUE 

_ WRITE  OUTPUT  TAPE  6*710 

WRITE  OUTPUT  TAPE  9,710 
71C  FCRHATIIOHR  •  ) 

N*  1 

DO  75C  W- 1 » I SL 
IF (IGNORE) 710 1,7 1C  1,7102 

..  7.1C  1  .  continue _ _ 

READ  INPUT  TAPES » 75 i^PI M  ),S(N),S(N+l),S(N+2)  ,5  IN *31 »  S I  N*4 ) , 
1  SIN  +  5  )  ,  S ( N+6 ) *  S  C  N ♦ 7 ) »  $  I N*  C ) »  S  I  N*9 ) 

7102  CONTINUE 
N*N-»LMAX 
KP=P(M) 

_  .  WRITE  OUTPUT  TAPE  6*75,2, M,KP 

WRITE  OUTPUT  TAPE  9,752,M»KP 


751 

FORMAT  IF1C.C, 1CF5 

.0) 

752 

FORMAT (15HS 

15,110 

75C 

CONTINUE 

WRITE  OUTPUT  TAPE 

6,753 

WRITE  OUTPUT  TAPE 

9,753 

'  753* 

FCRMATIllHS 

•  ) 

K  =  LMAX  *  I SL 

DC  754  N= 1 , K 

KS*XFIXF(S(N) ) 

WRITE  OUTPUT  TAPE 

6 , 755 , N , K S 

WRITE  OUTPUT  TAPE 

9,755, N,KS 

"7*55  ~* 

FORMAT  1 1 5HD 

15,110 

754 

CONTINUE 

WRITE  OUTPUT  TAPE 

6,756 

WRITE  OUTPUT  TAPE 

9,756 

756 

FORMAT (11HD 

•  ) 

T4*T44  1 . 

II  =  II*  ~ 

DO  759  I =1 , I SL 
DC  759  J=1,ISL 
DC  759  L*1,LMAX 
OOP! I, J,L)*C. 

_  __  15 9  .  TERM  ( I ,J)*0. 

. : . — IGNORE  ■  1  -- - - - - - 

GC  TO  1 
714  CONTINUE 

WRITE  OUTPUT  TAPE  6,757 
WRITE  OUTPUT  TAPE  9,757 

_  _  757_  FORMAT (17H  _ENO_CF  _PR0_8LEM  I 

WRITE  OUTPUT  ~TAPE  *6,758  ** 

WRITE  OUTPUT  TAPE  9,758 
756  FORMAT  1 1HC ) 

ENO  FILE  9 
REWINC  9 

i _ SECT  I  CN  I  -  CALCULATION  OF  MOVEMENT,  TI WES 

_ _ WRJ  TE  OUTPUT  _TAP£  4,U3 

713  FORMAT  ( 18 H  PROBLEM  COMPLETED  ) 

_ _ CALI  EXIT  . _ _ ...  .  . 
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